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An  efficient  technique  for  mapping  thermal  boundary  conditions  is  described  and 
demonstrated.  The  technique  is  based  on  a  piece-wise  polynomial  approximation  where  the 
Laplacian  derivatives  in  space  are  constrained  using  the  heat  equation.  Measured  values  for 
the  Laplacian  are  obtained  from  temperature  rate  measurements  from  sensors  embedded 
within  a  body.  The  technique  has  been  implemented  in  a  digital  signal  processor  and  is  able 
to  provide  real-time  data  on  thermal  boundary  conditions  over  a  surface.  The  technique  is 
adaptable  to  complex  geometry.  In  this  paper  the  technique  will  be  applied  to  a  study  of  the 
injector-wall  interactions  in  a  laboratory  scale  liquid  rocket  engine5. 


Nomenclature 


at 

= 

a  constant  in  a  polynomial,  also 

d 

= 

vector  of  data 

G 

= 

a  matrix  of  polynomial  terms 

P 

= 

a  vector  of  polynomials 

q 

= 

heat  flux 

T 

= 

temperature 

t 

= 

time 

w 

= 

plate  thickness 

x,y,z 

= 

dimensional  coordinates 

o 

= 

a  complex  constant 

CJ 

= 

angular  frequency 

* 

= 

complex  conjugate 

~ 

= 

frequency  domain  variable 

T 

= 

transpose 

-1 

= 

inverse 

I.  Introduction 

Local  concentrations  of  heat  flux  and  temperature  are  of  interest  in  aerospace  applications  including  injector- wall 
interactions  in  liquid  rocket  engines,  shock-shock  interactions  in  external  flow  over  hypersonic  vehicles, 
internal  and  external  flow  over  blades  and  vanes  in  gas-turbine  engines,  and  deposition  of  directed  energy.  Accurate 
characterization  of  the  distributions  may  be  required  for  validation  of  a  model,  optimization  of  a  design,  or  control 
of  an  operational  system.  There  are  two  basic  methods  in  use  today.  In  the  calorimetry  approach  a  heat  exchanger  is 
incorporated  into  the  wall  of  the  device  and  the  temperature  rise  of  the  working  fluid  is  used  to  determine  the  local 
heat  flux.  This  method  is  limited  to  applications  where  the  heat  flux  is  slowly  varying  in  time  and  changes  in  one 
direction  only  along  the  surface.  Various  assumptions  and  approximations  are  required  to  infer  the  surface 
temperature;  although,  in  some  cases,  temperature  sensors  can  be  incorporated  to  give  direct  measurements.  To 
obtain  high  spatial  resolution  with  a  calorimeter,  specialized  manufacturing  processes  are  required.  The  other  major 
class  of  methods  is  appropriate  for  transient  situations  which  commonly  occur  in  hypersonic  wind  tunnel  test  articles 
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and  laboratory  scale  combustion  devices  where  the  test  articles  are  designed  as  heat-sink  devices.  Surface  junction 
thermocouples  have  been  developed  which  match  the  thermal  characteristics  of  the  test  articles  and  can  be  treated  as 
non-intrusive  to  the  flow  of  heat.  From  the  temperature  history  a  heat  flux  history  can  be  constructed.  In  some 
situations  the  temperature  sensors  must  be  embedded  in  the  wall  at  some  distance  from  the  surface.  In  this  situation, 
inverse  methods  must  be  used  to  reconstruct  the  surface  heat  flux  and  temperature.  The  mathematical  methods  for 
performing  the  analysis  for  surface  junction  and  embedded  sensors  have  been  extensively  studied  and  are  highly 
developed;  however  the  most  widely  used  methods  assume  that  the  flow  of  heat  is  one-dimensional.  The  potential 
errors  in  this  approach  can  be  readily  appreciated  when  one  considers  the  situation  that  occurs  when  heat  is 
deposited  by  a  concentrated  source  such  as  a  laser.  When  energy  is  deposited  adjacent  to  the  sensor,  heat  diffuses  in 
the  plane  of  the  surface  causing  a  temperature  rise  at  the  location  of  the  sensor  which  is  recorded  incorrectly  as  heat 
flux  at  that  point.  Due  to  a  plethora  of  such  applications,  research  in  this  area  has  shifted  to  the  development  of 
solution  techniques  for  spatially  multi-dimensional  and  transient  problems. 

The  mathematical  characteristics  of  the  inverse  heat  conduction  problem  are  well  known.  Changes  in 
boundary  conditions  cause  attenuated  and  time-lagged  responses  within  a  body.  In  the  presence  of  noise  and 
measurement  error,  the  problem  becomes  ill-posed  and  solutions  can  be  biased  and  unstable.  Successful  inverse 
techniques  have  been  developed  which  rely  on  the  use  of  information  from  future  time  steps  to  ensure  that  changes 
in  the  boundary  conditions  at  the  current  time  step  have  had  time  to  reach  the  sensor  locations.  Stabilization  has  also 
been  accomplished  through  Tikhonov  regularization.  Other  techniques  utilize  low-order  basis  functions  to  dampen 
instabilities.  Most  techniques  rely  on  an  iterative  solution  wherein  an  initial  estimate  for  the  boundary  condition  is 
made,  followed  by  a  forward  calculation  to  establish  the  effect  at  the  sensor  locations,  and  then  a  refinement  in  the 
boundary  condition  with  repetition  until  a  stopping  criteria  is  satisfied.  The  need  to  solve  the  forward  problem 
multiple  times  at  each  time  step  in  a  transient  analysis  results  in  lengthy  run  times  and  must  be  performed  in  post¬ 
processing.  In  addition  to  the  destabilizing  effects  of  noise  and  measurement  error  and  the  time  required  to  complete 
a  calculation,  practical  methods  must  also  contend  with  the  temperature  dependence  of  thermal  properties,  complex 
geometries,  the  need  to  accurately  model  the  effects  of  temperature  sensor  installations,  have  the  ability  to  handle 
transient  and  steady  state  applications,  and  be  intelligible  to  those  responsible  for  applying  the  technique. 


II.  Theory 

In  this  work  we  examine  the  use  of  the  finite  element  method.  Finite  element  methods  are  highly  developed 
for  forward  problems.  They  are  used  because  the  basic  approach  of  subdividing  the  domain  into  cells  and 
approximating  the  solution  over  each  cell  with  a  polynomial  is  readily  adaptable  to  complex  geometries  and  is 
mathematically  tractable.  However,  there  are  some  distinct  differences  between  forward  and  inverse  problems  that 
result  in  some  changes  in  the  basic  strategy.  In  solving  forward  problems  it  is  possible  to  generate  meshes  that  are 
limited  in  size  only  by  the  memory  capacity  of  the  computer.  It  is  possible  to  finely  discretize  regions  where  the 
function  behavior  is  complex  such  that,  within  each  individual  cell,  the  function  can  be  approximated  adequately 
with  a  simple  linear  basis  function.  In  inverse  problems,  as  developed  here,  the  grid  corresponds  to  a  network  of 
sensors.  In  practical  applications,  grids  with  tens  or  hundreds  of  nodes  are  possible,  but  probably  not  thousands,  and 
certainly  not  the  millions  that  are  now  common  in  forward  problems.  This  provides  a  strong  impetus  for  developing 
schemes  that  use  high  order  basis  functions.  Coy  showed  how  basis  functions  of  arbitrary  order  can  be  developed  in 
one  spatial  dimension  using  two  nodes  by  constraining  the  approximating  polynomials  using  the  even  numbered 
spatial  derivatives  from  the  heat  equation.  Adding  second  derivatives  improved  the  frequency  response  by  three 
orders  of  magnitude.  The  model  could  be  formulated  as  a  digital  filter  and  calculated  in  real  time.  In  this  paper  we 
extend  the  approach  to  three  dimensions. 

To  fix  ideas  we  consider  a  set  of  sensors  embedded  in  a  solid  material  with  smoothly  varying  thermal 
transport  properties.  We  are  primarily  interested  in  constructing  approximations  for  the  temperature  and  heat  flux 
distribution  at  the  surface.  The  temperature  distribution  at  any  point  in  time  within  each  finite  element  is 
approximated  by  a  polynomial. 

TCxvy, z..t)  &  aD  -F  aLx  -F  a^y  -F  a2z  -F  a±xy  -F  a5xz  -F  a6yz  +-  —  (j) 

Or  simply  as, 


T(jcf  t)  ftf  PJ  a 
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(2) 


Where, 


(3) 


P1  =  [1  x  ;y  z  xy  xz  yz  .«.] 


<3-  —  GZ- Q-j  ■flg  hi  ] 


(4) 


In  many  applications  the  placement  of  the  sensors  and  the  choice  of  approximating  polynomial  can  be 
optimized  based  on  knowledge  of  the  physical  problem  at  hand.  If  it  is  known  that  gradients  in  one  particular 
direction  are  much  greater  than  other  directions,  the  sensor  spacing  can  be  adjusted  and  high  order  terms  in  only  the 
direction  of  the  high  gradient  need  to  be  included.  If  we  assume, 

&2  T  & 2  T  & 2  T 

- » - £f - 

dz*  '  3v2  (5) 

The  heat  equation  becomes, 


dT  _  dz  T 

&t  G  dz2  (6) 

The  time  derivatives  can  be  obtained  from  the  temperature  records  at  the  measurement  points.  The  set  of  equations 
to  be  solved  can  be  written, 


Or  simply  as, 


i’VJ 

dz2 

PP'lXn) 


■  Tt  ' 

■ay 

Tr. 

■am- 

Vai. 

-fn/aa- 

dz2 


Ga  =  d 


(7) 

(8) 


If  the  number  of  equations  exceeds  the  number  of  terms  in  the  polynomial,  we  can  solve  as  a  least-squares  problem. 


a  =  (GTGj~1GTd 


(9) 


Note  that  (CT  G)  1  Gr  depends  only  on  the  location  of  the  sensors. 

The  performance  of  the  method  can  be  evaluated  by  comparing  with  analytical  solutions.  Cole  has 
developed  Green’s  function  solutions  for  steady-periodic  heating  that  are  suitable  for  this  purpose.  Consider  a 
uniform  slab  of  depth  W  and  infinite  width  with  a  uniform,  single-frequency,  steady-periodic  heat  flux  on  the 
surface  at  y=W,  -a<x<a. 
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The  temperature  inside  the  slab  can  be  expressed  as, 

T(x,  y,  fr)  =  Real  [f  Of,  y,  &>}  eJi)f  ] 


(10) 


,  ct?= 
T{x,y,bi]  =  — 


1  ,«  c-<*l*-*rl  2  frm  ^  r< 

—  I  - dx'  +■  /  (— l)n — CDS  ( — yl  I 

W  )_.  2 aaD  LjT  IV  MV y}  J_ 

r /?i7F 

^  =  [^J  +^"/° 


®  0-<=iik-'ri 
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■dx' 
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9^*  sinhC.Gu.J 


flflr 


,x  <  —a 


_£  2  E  tJp 


■  dx'  =  { 


,  _  -a*Ca+jj  _ 


2a  a  £ 


-a  <  x  <  a 


9  ***  sinhCaUp] 


-,x  >  □. 


We  will  use  the  analytical  function  to  obtain  inputs  to  the  finite  element  model,  and  again  to  compare  the 
amplitude  and  phase  of  the  analytical  function  with  predictions  of  the  finite  element  model  at  points  distant  from  the 
sensor  locations,  most  importantly  at  the  surface.  The  model  coefficients  are  real  constants  that  multiply  the  data 
vector  so  we  can  write, 


a  =  (CtC)-1Gt  3 

•  rL  ■ 

f 

a  =  *, 

The  complex  temperature  predicted  by  the  finite  element  model  can  be  written, 


(11) 


(12) 


f(*)=PTS 


(13) 
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The  complex  heat  flux  in  the  z-direction  is, 


gp'1 

qjr£UJ  ftf  —k  —  a 
oz 

The  gain  and  phase-lag  of  the  finite  element  model  relative  to  the  analytical  solutions  are  the  following, 

r{X,  frj)  -  j-.-p-  _ 

p  (*..&*)  =  atan  Umag  (TFE')/ Real  —  Gtan[lmag(f')/R&al(f')] 


'S, 

)  w 


*iL, 


< 


I  1  111! 

•  •  1 
•  •  i 

i 


(14) 


(15) 

(16) 


At  locations  x  =  ±a  there  are  gradients  of  temperature  in  the  x-direction  that  we  would  like  to  capture.  We 
take  as  a  length  scale  the  distance  from  the  heated  surface  to  the  first  sensor,  ff£.  Thus  we  have, 
=  1,  fi'3+  =  S2/S v  a+  =  a/Slfl+  —  l/SlfW+  =  WfSL,x+  =  x/6Lfz+  =z/St.  Temperature  and  frequency 
are  non-dimensionalized  as  follows:  f+  =  fk/qQa,  and  w+  —  w&l/a.  We  consider  a  case  where  o-+  and  W+  are 
large  and  look  at  the  effect  of  frequency,  ^+,  on  the  surface  temperature,  T+f  amplitude  and  phase  for  several  values 
of  the  sensor  spacing,  f+. 


III.  Application 

In  the  final  version  of  this  paper  we  will  include  information  on  how  we  have  implemented  the  technique 
using  a  digital  signal  processor  to  obtain  real-time,  video-like  imaging  of  thermal  boundary  conditions.  We  will  also 
show  an  example  application  in  the  measurement  of  thermal  boundary  conditions  in  a  laboratory  scale  rocket, 
showing  patterns  of  heat  flux  produced  by  a  multi-element  injector,  and  how  these  patterns  change  as  we 
deliberately  introduce  non-uniformities  into  the  flows  from  the  injector  elements. 
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